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Rotating BHs as Central Engines of Long GRBs: 

Faster is Better 



O 

(N 

£ 

1—5 

00 
(N 

43 

6 

CO 

03 



> 

On 



13 



Shigchiro Nagataki 

Yukawa Institute for Theoretical Physics, Kyoto University, Oiwake-cho Kitashirakawa Sakyo-ku, Kyoto 606-8502, Japan 

nagataki@yukawa. kyoto-u. ac.jp 

(Received ; accepted ) 
Abstract 

We performed simulations of collapsars with different Kerr parameters a =0, 0.5, 0.9, 0.95. It is shown 
that a more rapidly rotating black hole is driving a more energetic jet. No jet is seen for the case of 
Schwartzschild black hole case, while the total energy of the jet is as large as 10 50 erg for a rapidly rotating 
Kerr black hole case (a = 0.95). In order to explain the high luminosity of a GRB, it is concluded that a 
rapidly rotating black hole is favored ('faster is better'). We also find in the case of a = 0.95 that (i) the 
stagnation region is clearly found in the jet region, (ii) the ordered poloidal field lines are seen in the jet, 
(iii) the jet region is surrounded by a 'Wall-like' structure that has a higher pressure than the jet region 
and contains strong vertical magnetic fields, and (iv) the jet is initiated by outgoing Poynting flux from 
the outer horizon of the black hole (Blandford-Znajck effect). The bulk Lorentz factor of the jet is still of 
the order of unity. However, energy density of electro-magnetic fields dominates the one of rest-mass in 
the jet. It can be expected that a relativistic jet will be seen if we perform a simulation for a longer time 
scale (of the order of 10 — 100 sec). 

Key words: gamma rays:bursts — black hole physics — relativity — supernovac: general — accretion, 
accretion disks 



1. Introduction 

It is still unknown how the central engine of Long 
Gamma-Ray Bursts (hereafter it is called as GRBs for 
simplicity) is working at the center of massive stars. In 
other words, we poorly know the engine that drives the 
most powerful explosion in the universe. Some of the su- 
pernovae that associate with GRBs were very energetic. 
Their energies were of the order of 10 52 ergs, which is 
~ 10 times larger than the ones of normal core-collapse 
supernovae. This fact strongly suggests that other engine 
from the one of normal core-collapse supernovac should 
be working at the center of the GRB's progenitors. 

There are some possible scenarios to drive a GRB jet. 
One of the most promising scenarios is the collapsar sce- 
nario (Woosley 1993). In this study, we investigate the 
collapsar scenario. In the collapsar scenario, a black hole 
is formed as a result of gravitational collapse. Also, rota- 
tion of the progenitor plays an important role. Due to the 
rotation, an accretion disk is formed around the equato- 
rial plane. After the formation of the accretion disk, there 
are several possible stories as below. 

One is that a jet-induced explosion along the rota- 
tion axis may occur due to the heating through neu- 
trino anti-neutrino pair annihilation that are emitted 
from the accretion disk (Woosley 1993; MacFadyen & 
Woosley 1999; Fryer & Meszaros 2000; Nagataki ct al. 
2003a; Nagataki et al. 2007; Sckiguchi and Shibata 
2007; Barkov and Komissarov 2010; Harikae et al. 2010). 
This is an interesting possibility, although neutrino heat- 
ing efficiency looks small (0.1 — 1%). A very careful treat- 



ment of neutrino transfer with conserved-scheme of hy- 
drodynamics will be necessary to prove that this effect is 
the key-process of GRB's engine. 

Another is that a jet is driven by extracting rotational 
energies of the accretion disk with a help of magnetic fields 
that pierce the disk (Blandford-Paync effect: Blandford 
and Payne 1982). This scenario is also investigated by 
several authors (Proga et al. 2003; Proga and Begelman 
2003; Mizuno et al. 2004a; Mizuno et al. 2004b; Proga 
2005; Fujimoto et al. 2006; Nagataki et al. 2007; Suwa 
et al. 2007; Harikae et al. 2009). This is also a promis- 
ing possibility. Further study will be necessary whether a 
relativistic jet will be launched by this mechanism. 

Recently, the effect of extraction of rotation energy from 
the black hole through outgoing Poynting flux (Blandford- 
Znajck effect: Blandford and Znajek 1977) is inves- 
tigated by using a General Relativistic Magneto-hydro 
Dynamics (GRMHD) code (Barkov and Komissarov 
2008; Komissarov and Barkov 2009a; Nagataki 2009; 
Barkov and Komissarov 2010). Energy extraction from 
a rotating black hole is a general relativistic effect, so 
GRMHD code is necessary to investigate the effect. In 
Barkov and Komissarov (2008), they showed that a jet is 
launched by Blandford-Znajek effect using a Kerr black 
hole with Kerr parameter (a=0.9) and a polytrope den- 
sity profile for a massive star model. In their succes- 
sive papers (Komissarov and Barkov (2009a), Barkov and 
Komissarov 2010), they also succeeded to launch a jet by 
Blandford-Zna jek effect (Kerr Parameter was chosen to be 
0.9 in Komissarov and Barkov 2009a, while 0.46 and 0.6 in 
Barkov and Komissarov 2010) using a polytrope density 
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profile for a massive star model. In Nagataki (2009), it 
was shown that a jet is launched by Blandford-Znajek ef- 
fect (Kerr parameter was chosen to be 0.5) using a realistic 
progenitor model developed by Wooslcy and Hegcr (2006). 
In Nagataki (2009), it was pointed out that Blandford- 
Payne effect may be also working. 

As stated above, in the previous papers, it was shown 
that a jet is successfully launched by Blandford-Znajek 
effect. However, there is no systematic study how the dy- 
namics depends on the Kerr parameter. It is true that 
Komissarov and Barkov (2010) presented results for dif- 
ferent Kerr parameters, but they also changed the ampli- 
tude of the initial magnetic fields. Also, no simulation has 
been reported for a = 0, that is, Schwartzschild black hole 
case. By performing a simulation for that case, we can 
clearly see how effectively the rotating black hole is work- 
ing to drive the jet. In this study, we present 4 simulations 
for the same initial condition with Nagataki (2009), but 
with different Kerr parameters a =0, 0.5, 0.9, 0.95. It is 
shown clearly that a more rapidly rotating black hole is 
driving a more energetic jet. Especially, in the case of the 
Schwartzschild case, no jet is found. That proves that the 
jet is really driven by the rotating black hole. 

In section 2, method of calculation is explained. In 
section 3, initial condition is shown. Results are given in 
section 4. Discussion is presented in section 5. Conclusion 
is stated in section 6. 

2. Method of Calculation 

In this study, we use the GRMHD code developed in 
Nagataki (2009). Thus we briefly explain the method of 
calculation here, and please see Nagataki (2009) for de- 
tails. 

In Nagataki (2009), we have developed a two- 
dimensional GRMHD code following Gammie et al. 
(2003) and Noble et al. (2006). We have adopted a 
conservative, shock-capturing scheme with Harten, Lax, 
and van Leer (HLL) flux term (Harten et al. 1983) and 
flux-interpolated constrained transport technique (Toth 
2000). We use a third-order Total Variation Diminishing 
(TVD) Runge-Kutta method for evolution in time, while 
monotonized central slope-limited linear interpolation 
method is used for second-order accuracy in space (van 
Leer 1977). 2D scheme (2-dimensional Newton- Raphson 
method) is usually adopted for transforming conserved 
variables to primitive variables (Noble et al. 2006). 

When we perform simulations of GRMHD, Modified 
Kerr-Schild coordinate is basically adopted with mass 
of the black hole (M) fixed where the Kerr-Schild ra- 
dius r is replaced by the logarismic radial coordinate 
x\ = lnr. When we show the result, the coordinates are 
transfered from Modified Kerr-Schild coordinates to Kerr- 
Schild ones. In the following, we use G = M = c = 1 unit. 
G is the gravitational constant, c is the speed of light, 
and M is the gravitational mass of the black hole at the 
center. 



3. Initial Condition 

The initial condition is also same with Nagataki (2009), 
but for different Kerr Parameters. Thus we explain it 
briefly here, and please see Nagataki (2009) for details. 

The calculated region corresponds to a quarter of the 
meridian plane under the assumption of axisymmetry and 
equatorial symmetry. The spherical mesh with 256 (r)x 
128(0) grid points is used for all the computations. The 
radial grid is nonuniform, extending from r =0.98r + to 
3xl0 4 (r+ = 1 + \/l — a 2 is the outer horizon, and 3xl0 4 
corresponds to 8.9x 10 9 cm in cgs units) with uniform grids 
in the Modified Kerr-Schild coordinate. 

We adopt the model 12TJ in Woosley and Heger (2006) 
(as for the progenitor evolution model, see e.g., Tutukov 
and Fedorova 2007). This model corresponds to a star 
that has 12Af Q initially with 1% of solar metallicity, and 
rotates rapidly and does not lose its angular momentum 
so much by adopting small mass loss rate. As a result, 
this star has a relatively large iron core of 1.82M Q , and 
rotates rapidly at the final stage. We assume that the 
central part of the star has collapsed and formed a black 
hole of 2M . Since M = 2M Q , r = 1 corresponds to 2.95 
x 10 5 cm. We also assume that the gravitational mass of 
the black hole is unchanged throughout the calculation. In 
Nagataki (2009), the Kerr parameter, a, was assumed to 
be 0.5, but in this study we perform simulations for a = 
0, 0.5, 0.9, 0.95, respectively (we name them as Model 
A,B,C,D, respectively). 

Since 1-D calculation is done for the model 12TJ, we 
can use the data directly for the physical quanta on the 
equatorial plane. As for the density, internal energy 
density, and radial velocity, we assume the structure of 
the star is spherically symmetric. We also set u 6 = 
initially. As for u^, we extrapolate its value such as 
u^(r,6») = ii^(r, tt/2) x sin6». 

We assume the vector potential oc max(p//9 max — 
0.2,0) sin 4 where /9 max is the peak density in the pro- 
genitor (after extracting the central part of the progenitor 
that has collapsed and formed a black hole). The field is 
normalized so that the minimum value of p gas /p mag = 10 2 
where p gas is the thermal pressure and p m ag is the mag- 
netic pressure. 

We use a simple equation of state p gas = (F — l)u where 
u is the internal energy density. We set T=4/3 so that the 
equation of state roughly represents radiation gas. 

As for the boundary condition in the radial direction, 
we adopt the outflow boundary condition for the inner 
and outer boundaries. As for the boundary condition in 
the zenith angle direction, axis of symmetry boundary 
condition is adopted for the rotation axis, while the re- 
flecting boundary condition is adopted for the equatorial 
plane. As for the magnetic fields, the equatorial symme- 
try boundary condition, in which the normal component 
is continuous and the tangential component is reflected, 
is adopted. 
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4. Results 

In figure 1, contours of rest mass density in logaris- 
mic scale for all models at the same time-slice t = 160000 
(that corresponds to 1.5760 sec) are shown. Cgs units 
are used for the rest mass density, while the length in the 
vertical/horizontal axes is written inG = M = c = l unit. 
r = 1 and 4000 corresponds to 2.95 xlO 5 cm and 1.18 
x 10 9 cm, respectively. These results are projected on the 
(r sin#, r cos #)-plane. Upper left panel shows the state 
of Model A (a = 0), upper right panel shows the one of 
Model B (a = 0.5), lower left panel shows the one of Model 
C (a = 0.9), and lower right panel shows the one of Model 
D (a = 0.95). It is clearly shown that the rotating black 
hole drives the jet (the Schwartzschild black hole cannot 
drive a jet (Model A), while a more rapidly rotating black 
hole is driving a stronger jet). 

From figure 2 to 6, we show the results for Model D at 
t = 160000 (that corresponds to 1.5760 sec). The length 
in the vertical/horizontal axes is written in G = M = c — 1 
unit (r = 1,20,100 correspond to 2.95 xlO 5 cm, 5.9 xlO 6 
cm, 2.95 xlO 7 cm, respectively). 

In figure 2, color-contours of rest mass density in log- 
arismic scale (cgs units) with velocity fields (arrows) are 
shown. It is noted that stagnation region where the ra- 
dial velocity becomes positive (outgoing) from negative 
(accreting) is seen around r = 15 in the jet. In contrast, 
it is shown later (in figure 5) that outgoing Poynting flux 
is positive from the outer horizon. 

In figure 3, color-contours of rest mass density in log- 
arismic scale (cgs units) with contours of the <fi compo- 
nent of the vector potential (A^) are shown. Level sur- 
faces coincide with poloidal magnetic field lines, and field 
line density corresponds to poloidal field strength. Upper 
panel shows the central region (20 x 20), while lower panel 
shows the wider region (100 x 100). In the upper panel, 
the ordered poloidal field lines are seen in the jet. In the 
lower panel, a 'Wall-like' structure that contains vertical 
magnetic fields is seen. 

In figure 4, contours of total pressure (sum of thermal 
and magnetic pressure) in logarismic scale. Cgs units are 
used for the pressure contours are shown. It is clearly 
seen that the 'Wall-like' structure has high pressure: it is 
higher than the one in the jet region. 

In figure 5, contours of outgoing Poynting flux (Fe in 
Eq.l) in logarismic scale. It is noted that the total energy 
flux, which is the integrated outgoing Poynting flux over 
the zenith angle, can be written as 

E = 2-k( d<V=^(-T£ Mi4 ) = 2tt / d0F E , (1) 
Jo Jo 

where g and T^ M t are determinant of the metric and 
(r,t) component of energy-momentum tensor of electro- 
magnetic fields. The unit of the contours is 10 50 erg s _1 
sr _1 . Outer horizon is seen at the center (r + = 1.312). 
It is clearly seen that outgoing Poynting flux is coming 
out from the outer horizon. This figure clearly shows that 
Blandford-Znajek process is working in this system. We 
could sec also a powerful outgoing Poynting flux toward 
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the accretion disk. Partially this may be working for the 
variability of the mass accretion rate seen in Nagataki 
(2009). 

In figure 6, contours of the electro-magnetic field energy 
flux per unit rest-mass flux are shown, which represent the 
bulk Lorentz factor of the invischid fluid element when 
all of the electro-magnetic field energy are converted into 
kinetic energy (Nagataki 2009). Even though the bulk 
Lorentz factor of the jet is still low (of the order of unity), 
the terminal bulk Lorentz factor can be relativistic. 

Finally, in figure 7, plots of the jet energy for all models 
at t = 160000 are shown. The definition of the jet energy 
is: 

/■OO r-Q 

E Jet = 2 X 2tt / dr d0y/=jj(Tt - pu°u ), (2) 

Jr + Jo 

where T\ is the (t,t) component of total energy- 
momentum tensor and integration is done only for the 
region where u r (radial component of 4-velocity of fluid) 
is positive. It is noted that the contribution of the rest 
mass energy is subtracted. Factor 2 is coming from the 
symmetry of the system with respect to the equatorial 
plane. Blue curve represents the jet energy within the 
opening angle 6> = 5°, while red curve represents the one 
within 6 — 10°. The unit of vertical axis is 10 48 erg. It is 
clearly seen that a more rapidly rotating black hole is driv- 
ing a stronger jet. The total energy of the jet for Model D 
(a = 0.95) is as large as 10 50 erg. For comparison, the ana- 
lytic curves of the BZ-flux formulation of Tchckhovskoy et 
al. (2010) for the case of monopole solution B=5 x 10 14 G 
are shown by the red-dashed curve (until forth power of 
= a/2rn where th is the horizon of the black hole). 
The analytic curve of the BZ-flux formulation of Tanabe 
and Nagataki (2008) for the case of B=5 x 10 14 G is shown 
by the black-dashed curve (until forth power of a). It is 
seen that the jet energy obtained in this study fits well 
with the previous analytical formulations. 

5. Discussion 

1. Whole Picture. 

In section 4, we could see that a stronger jet is driven by 
a more rapidly rotating black hole. It is clear that the 
rotating black hole is driving the jet. In fact, no jet is 
seen for the case of Model A (a = 0). On the other hand, 
the total energy of the jet for Model D (a = 0.95) is as 
large as 10 50 erg. We showed some figures for Model D 
to understand the dynamics more clearly. We have found 
that (i) the stagnation region is seen around r = 15 in the 
jet region, (ii) the ordered poloidal field lines are seen in 
the jet (that is consistent with other previous study on the 
jet formation from a torus (e.g. McKinney & Gammic 
2004)), (iii) the jet region is surrounded by a 'Wall-like' 
structure that has a higher pressure than the jet region 
and contains strong vertical magnetic fields, and (iv) the 
jet is initiated by outgoing Poynting flux from the outer 
horizon, effect (Blandford-Znajek effect). 

2. Wall-Like Structure. 

We found the 'Wall-like' structure around the jet. It is 
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considered that the structure will depend on the initial 
angular momentum and density profile of the progenitor 
star and initial configuration and its strength of magnetic 
fields. We should investigate this feature by performing 
some simulations systematically as a next step. However, 
this Wall-like structure seems not to work so much for the 
formation and collimation of the jet, because the colli- 
mated jet has been launched already before the formation 
of the Wall-like strucuture. 

3. Jet Energy. 

It has been shown that the jet energy can be as large as 
10 50 erg in 1.5 sec for Model D. Thus we expect that the 
energy of the jet can be ~ 10 51 erg in 10 sec, and we can 
say that Blandford-Znajek mechanism is very promising to 
drive a GRB jet. On the other hand, for a slowly-rotating 
black hole case, the jet power looks too low, because the 
luminosity of a GRB jet (10 51 erg s _1 ) is hard to explain. 
Thus it is concluded that a rapidly rotating black hole is 
favored ('faster is better'). This result is consistent with 
some previous test calculations (Tanabe and Nagataki 
2008; Tchekhovskoy et al. 2010). Due to the limit of 
computational resources, we had to stop the simulation 
at t = 160000 in this study. We are going to present a 
simulation that should last more than 10 sec in the near 
future. 

4. Lorentz Factor. 

It is found that the bulk Lorentz factor of the jet can- 
not be so high in this study (of the order of unity), 
although the electro-magnetic field energy flux is much 
greater than rest-mass flux in the jet, which was also seen 
in Nagataki (2009). We believe this jet has a good prop- 
erty, because the jet speed will become relativistic if the 
electro-magnetic field energy is converted to kinetic en- 
ergy through the propagation. Also, we can be optimistic 
because the Wall-like structure was found in our simu- 
lations. It is pointed out that such a Wall is necessary 
to drive a relativistic jet in some papers (Komissarov et 
al. 2009b; Tchekhovskoy et al. 2010). Moreover, there 
have been many papers on the propagation of relativistic 
jet from a massive star (Aloy et al. 2000; Aloy et al. 
2002; Zhang et al. 2003; Zhang et al. 2004; Mizuta 
et al. 2006; Mizuta and Aloy 2009; Morsony et al. 
2007; Wang et al. 2008; Morsony et al. 2010; Mizuta et 
al. 2010; Nagakura et al. 2010) where the jet is driven by 
injecting energy from the inner boundary. Some of their 
results suggest that the bulk Lorentz factor of the jet will 
be increasing with time (e.g. Mizuta et al. 2010) because 
the material of the progenitor star is blown off and the jet 
region becomes 'cleaner' with time. Thus we expect that 
a relativistic jet will be seen if we perform a simulation 
for a longer time scale (of the order of 10 — 100 sec. see 
also Suwa and Ioka 2010). 

5. Blandford-Payne Effect. 

In this study, we showed that Blandford-Znajek effect is 
working. On the other hand, it was found that no jet 
was driven with a help of rotation of the accretion disk 
(Blandford-Payne effect) in Model A (a = 0). Of course 
there is a possibility that the disk will drive a jet later. 
Also, the interaction between the rotating black hole and 



the accretion disk may be important. That is, if the 
rotating black hole is connected with the accretion disk 
by magnetic fields and rotation power is conveyed to the 
disk, Blandford-Payne effect will become more effective. 
Further study will be necessary to conclude the necessity 
of Blandford-Payne effect to drive a GRB jet. 

6. Microphysics. 

In this study, microphysics such as neutrino physics 
and nuclear reaction is not included. The disk struc- 
ture is deformed by neutrino cooling process (Popham 
et al. 1999; Nagataki and Kohri 2002; Nagataki et al. 
2003a; Nagataki et al. 2007; Lopez-Camara ct al. 
2010; Linder et al. 2010; Taylor et al. 2010; Sekiguchi 
and Shibata 2010). Nucleosynthesis study provides us 
with an opportunity to compare the simuations with ob- 
servations of supernova ejecta and supernova remnants 
(Nagataki et al. 1997; Nagataki ct al. 1998a; Nagataki et 
al. 1998b; Nagataki 2000; Maeda et al. 2002; Nagataki 
et al. 2003b; Takiwaki et al. 2004; Nagataki ct al. 
2006; Maeda et al. 2007; Tominaga 2009; Ramirez-Ruiz 
andMacFadyen 2010; Milosavljevic et al. 2010). We are 
planning to include these microphysics as a next step. 

7. 3D Simulation. 

We showed our results of two-dimensional (2D) simula- 
tions in this study. It will be important to compare the 
2D simulations with 3D ones. In 3D simulations, hy- 
drodynamic instability in the azimuthal direction will be 
seen (Nagakura & Yamada 2008; Nagakura & Yamada 
2009; Mckinney & Blandford 2009). 

6. Conclusion 

We have performed 4 simulations of collapsars with dif- 
ferent Kerr parameters a =0, 0.5, 0.9, 0.95. It is clearly 
shown that a more rapidly rotating black hole is driving 
a more energetic jet. No jet is seen for the case of Model 
A (a = 0), while the total energy of the jet for Model D 
(a = 0.95) is as large as 10 50 erg. We expect that the en- 
ergy of the jet can be ~ 10 51 erg in 10 sec for Model D. 
In order to explain the high luminosity of a GRB, it is 
concluded that 'faster is better'. 

We have found in Model D that (i) the stagnation region 
is seen around r = 15 in the jet region, (ii) the ordered 
poloidal field lines are seen in the jet, (iii) the jet region is 
surrounded by the 'Wall-like' structure that has a higher 
pressure than the jet region and contains strong vertical 
magnetic fields, and (iv) the jet is initiated by outgoing 
Poynting flux from the outer horizon (Blandford-Znajek 
effect). 

As for the bulk Lorentz factor, it is still not so high in 
this study (of the order of unity) , although energy density 
of electro-magnetic field dominates the one of rest-mass 
in the jet. We expect that a relativistic jet will be seen 
if we perform a simulation for a longer time scale (of the 
order of 10 — 100 sec). 
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Fig. 1. Contours of rest mass density in logarismic scale for 
all models at the same time-slice t = 160000 (that corresponds 
to 1.5760 sec). Cgs units are used for the rest mass density, 
while the length in the vertical/horizontal axes is written in 
G = M = c = 1 unit, r = 1 and 4000 corresponds to 2.95 
xlO 5 cm and 1.18 xlO 9 cm, respectively. These results are 
projected on the (r sin0, r cos 0)-plane. Upper left panel 
shows the state of Model A (a = 0), upper right panel shows 
the one of Model B (a = 0.5), lower left panel shows the one 
of Model C (a = 0.9), and lower right panel shows the one of 
Model D (a = 0.95). 
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Fig. 3. Same with Figure 1, but for Model D (a = 0.95) with 
line-contours of the <f> component of the vector potential (A$). 
Upper panel shows the central region (20x20inG = M = c= l 
unit), while lower panel shows the wider region (100 X 100). 
r = 100 corresponds to 2.95 xlO 7 cm. 
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Fig. 4. Contours of total pressure (sum of thermal and mag- 
netic pressure) in logarismic scale at t = 160000 (that corre- 
sponds to 1.5760 sec) for Model D (a = 0.95). Cgs units are 
used for the pressure contours, while the length in the verti- 
cal/horizontal axes is written inG = M = c= l unit, r = 100 
corresponds to 2.95 XlO 7 cm. 




Fig. 6. Contours of the energy flux per unit rest-mass flux 
for Model D (a = 0.95) at t = 160000. The contours represent 
the bulk Lorentz factor of the invischid fluid element when all 
of the internal and magnetic energy are converted into kinetic 
energy. 




Fig. 5. Contours of outgoing Poynting flux in logarismic 
scale at t = 160000 (that corresponds to 1.5760 sec) for Model 
D (a = 0.95). The unit of the contours is 10 50 erg s — 1 sr -1 
(see text for the definition in detail), while the length in the 
vertical /horizontal axes is written inG = M = c= l unit, r = 1 



corresponds to 2.95 
center (r+ = 1.312). 



xlO cm. Outer horizon is seen at the 
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Fig. 7. Plots of the jet energy (see text for the definition in 
detail) for all models at t = 160000 (that corresponds to 1.5760 
sec). The unit of vertical axis is 10 48 erg. Blue curve repre- 
sents the jet energy within the opening angle 9 = 5° , while red 
curve represents the one within 8 = 10° . The analytic curves 
of the BZ-flux formulation of Tchekhovskoy et al. (2010) for 
the case of monopole solution with B=5 X 10 14 G are shown by 
the red-dashed curve (forth power of Qn ) . The analytic curve 
of the BZ-flux formulation of Tanabe and Nagataki (2008) for 
the case of monopole solution with B=5 X 10 14 G is shown by 
the black-dashed curve (until forth power of a). 



